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We exploit the renormalization-group approach to estabhsh the exact infrared behavior of an 
interacting Bose system at zero temperature. The local-gauge symmetry in the broken-symmetry 
phase is implemented through the associated Ward identities, which reduce the number of indepen- 
dent running couplings to a single one. For this coupling the e-expansion can be controlled to all 
orders in e (= 3 — d). For spatial dimensions 1 < d < 3 the Bogoliubov fixed point is unstable towards 
a different fixed point characterized by the divergence of the longitudinal correlation function. The 
Bogoliubov linear spectrum, however, is found to be independent from the critical behavior of this 
correlation function, being exactly constrained by Ward identities. The new fixed point properly 
gives a finite value of the coupling among transverse fluctuations, but due to virtual intermediate 
longitudinal fluctuations the effective coupling affecting the transverse correlation function flows to 
zero. As a result, no transverse anomalous dimension is present. This treatment allows us to recover 
known results for the quantum Bose gas in the context of a unifying framework and also to reveal 
the non-trivial skeleton structure of its perturbation theory. 



I. INTRODUCTION 



Perturbation theory for neutral zero temperature in- 
teracting bosons in the superfluid phase has been devel- 
oped in the fifties"'^'^ following the mean-field result of 
Bogoliubov.'^ Notwithstanding the success of its imme- 
diate implementation in the low density limit ^ and in 
proving the coincidence of the quasiparticle and sound 
modes, the unresolved problem of infrared (IR) diver- 
gences due to the Goldstone mode led to surprising re- 
sults. As a matter of fact, much later A. Nepomnyashchii 
and Y. Nepomnyashchii^ (from here on referred as NN) 
demonstrated that the fundamental quantity I]i2(fc), the 
anomalous self-energy, must actually vanish for vanish- 
ing momentum {k — > 0) and d < 3. In the Bogoliubov 
approximation the existence of a linear spectrum relies 
on the presence of Ei2(fc) ^ 0. Since NN found that 
the linear spectrum is preserved for A; — > 0, their result 
implies that the structure of the quadratic part of the ac- 
tion must be actually much more complicated than the 
simple Bogoliubov form. 

The NN results were obtained by studying directly 
the skeleton structure of the diagrammatic perturbation 
theory (PT). To determine the asymptotic behavior of 
the correlation functions for A: ^ NN had to resort 
to a self-consistent analysis of the skeleton perturba- 
tion theory, finding relations among set of diverging dia- 
grams. This method, besides being not fully transparent, 
might not be readily controlled. The need of develop- 
ing a PT free from IR divergences was actually realized 
by Popov^ before NN. Popov development was achieved 
within a functional-integral approach, by introducing a 
phase-amplitude representation for the bosonic variables 
in the IR region for |k| smaller than a cutoff kg- In this 



way at all intermediate steps the calculation turned out 
to be free from IR divergences. However, elimination of 
the arbitrary parameter kg remains nontrivial in Popov 
approach, as it requires full control of the IR region up to 
momenta of the order of ko where the phase-amplitude 
description becomes less and less accurate since at large 
momenta it is the particle representation to be appropri- 
ate. In addition, the phase-amplitude description some- 
what obscures the particle viewpoint even at small mo- 
menta. As a matter of fact, when A: it is not easy 
to follow with the phase-amplitude description the non- 
trivial evolution of the elementary excitations from free 
particles to the sound mode. Probably, it was for this rea- 
son that Popov recognized in his own language the NN 
result about the vanishing of the anomalous self-energy 
only at a later stage. 

The vanishing of the anomalous self-energy is not a 
mere mathematical result but it has a definite physi- 
cal origin. In this particular context of broken gauge 
symmetry, this finding reflects the general picture given 
by Patasinskij and Pokrovskij,^^ whereby the Goldstone- 
mode singularity of the transverse correlation function 
drives a divergence also in the longitudinal correlation 
function for a continuous broken symmetry. An addi- 
tional interesting issue is to understand how, technically, 
the IR divergences do not lead to a anomalous dimension 
for the quadratic propagator, as it happens instead for 
critical phenomena. 

These circumstances lead to the necessity of provid- 
ing a unifying and fully controlled treatment of IR di- 
vergences. In this respect, the renormalization group 
(RG) approach appears to be the natural and reliable 
tool to determine the IR behavior of the system in the 
presence of IR divergences. In this paper, we exploit 
the RG approach to obtain the exact IR behavior of the 
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vertex functions (and thus of the correlation functions) 
for a neutral system of interacting bosons in the broken- 
symmetry phase at zero temperature. To this end, we 
will make extensive use of the Ward identities (WI) as- 
sociated with gauge invariance, which pose strong con- 
straints on the RG equations and enable us to obtain 
the desired solution to all orders of the e-expansion (here 
e = 3 — d). Although Popov method might at a first look 
be preferable to deal with the broken-symmetry phase 
since it deals at the outset with variables that are man- 
ifestly gauge invariant, it appears useful to complement 
the phase-amplitude approach by studying the IR behav- 
ior of the system via a more standard particle represen- 
tation. Our treatment allows us to follow the evolution 
of the propagator continuously, from the particle region 
to the hydro dynamical one. A short account of this work 
has been given in Ref. 12. 

Usually, the appearance of IR divergences is related to 
a second-order phase transition, when a competition be- 
tween two different phases with equal free energy leads 
to divergent fluctuations at the critical temperature. In 
the present case, although the system is in a stable phase 
and no phase transition occurs, there exists a competi- 
tion among degenerate states which are associated with 
different values of the macroscopic phase of the order pa- 
rameter and are thus physically equivalent. This leads 
to the divergence of the longitudinal correlation function 
mentioned above, while all correlation functions obtained 
as averages of local gauge-invariant operators (like the 
density-density correlation function) are expected to be 
free from IR divergence. A similar behavior is found for 
fermionic systems, where the perturbation theory for the 
stable phase of the Luttinger liquid presents infrared di- 
vergences that have to disappear when the symmetry is 
correctly enforced through Ward identities. 

Perturbation theory for the broken symmetry phase in 
the particle representation considers averages of opera- 
tors which are not local gauge invariant, and is therefore 
plagued by IR divergences at intermediate steps of the 
calculation. It is clear that these divergences cannot be 
independent from each other as they have to cancel out 
when calculating averages of the above-mentioned local 
gauge-invariant operators. In this respect. Ward identi- 
ties provide explicit connections among the divergences. 
In turn, this makes all local gauge-invariant susceptibil- 
ities (like the compressibilities) finite and stabilizes the 
system accordingly with respect to phase fluctuations. 

By our RG approach we will be able to prove that 
the three available gauge-invariant susceptibilities of the 
system (namely, the condensate susceptibility, the ordi- 
nary susceptibility related to the sound velocity, and the 
total density) are indeed finite, being invariants of the 
RG flow. This will be done by establishing explicit con- 
nections among the running couplings of the theory. In 
this way, we will be left with only one independent run- 
ning coupling, whose infrared behavior will then be es- 
tablished to all orders in e (= 3 — d) for d > 1. The 
flxed point for < 3 is characterized by a finite inter- 



action coupling among the transverse fluctuating fields. 
This at first sight suggests that the theory is interacting 
at the fixed point like the standard 0^ theory. We find, 
in fact, that this is not true due to the presence of the 
longitudinal field. One can define an effective interac- 
tion among transverse fiuctuations entering the relevant 
subset of diagrams responsible of the appearance of the 
anomalous dimension in 0^ theory. This effective interac- 
tion is flowing to zero at the fixed point, thus preserving 
the behavior for the transverse propagator. The 

occurrence of the sound spectrum will turn out to be in- 
dependent from the scaling behavior of the independent 
running couplings, as it depends only on the underlying 
gauge symmetry. We will also recover the leading IR be- 
havior of the single-particle Green's functions and of the 
(two-particle) response functions obtained previously by 
Gavorct and Nozieres'' and by NN**. 

RG treatments of the zero-temperature interacting- 
boson problem have been previously given by Weich- 
man^''' and by Benfatto^*'. Weichman found the correct 
exponent of the longitudinal susceptibility by performing 
a one-loop calculation for an intrinsically non-divergent 
quantity (the free energy) and then used a simple scaling 
argument to determine the singular behavior of the longi- 
tudinal single-particle Green's function. This scaling ar- 
gument can, in general, be justified by the RG approach, 
as it will be shown below by our treatment. Benfatto, 
on the other hand, used a Wilson- like RG approach in 
a rigorous fashion to determine the scaling behavior for 
a number of running couplings. However, his treatment 
was limited to d = 3 and did not take advantage of a 
systematic implementation of WI, which in turn would 
guarantee the linear behavior of the spectrum irrespec- 
tive of the IR behavior of the longitudinal correlation 
function. The fact that the theory is asymptotically free 
in = 3, as found by Benfatto, is thus not essential 
to stabilize the superfluid phase by preserving the linear 
spectrum. 

The question naturally arises whether the momentum 
region, where phase fiuctuations (leading to IR diver- 
gences) dominate, has physical relevance apart from es- 
tablishing the superfiuid behavior. To this end, a general- 
ized Ginzburg criterion can be introduced for the momen- 
tum variable to determine, in particular, whether the IR 
region will or will not merge into the Bogoliubov region, 
where the spectrum is still linear but the longitudinal 
correlation function is approximately constant. For the 
low-density Bose gas one can explicitly show^° that the 
Ginzburg region extends up to a characteristic momen- 
tum that vanishes exponentially (in the gas parameter) 
with respect to the extension of the linear Bogoliubov re- 
gion. In the more general case, it is not a priori possible 
to assess whether the Bogoliubov region might be washed 
out by the growing of the Ginzburg region. In this case 
a full RG approach is required, as shown in the present 
paper. 

The plan of the paper is as follows. In Section II we 
recall how IR divergences manifest themselves in FT by 
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considering the one-loop correction to the Bogoliubov ap- 
proximation, and then set up the RG treatment for the 
problem at hand. In Section III we derive and analyze 
the WI relevant to our purposes. In Section IV we de- 
termine the exact IR behavior of the vertex functions in 
the context of the e expansion. Section V gives our con- 
clusions. For completeness, additional technical points 
are discussed in the Appendices. Specifically, Appendix 
A deals with the definition of gauge invariancc in terms 
of the couplings introduced in this paper, Appendix B 
shows the one-loop evaluations of relevant vertex func- 
tions and the validity of the WI, Appendix C gives the 
one-loop equations for the running couplings entering the 
propagators. 



II. IR DIVERGENCES AND SETUP OF THE RG 
TREATMENT 

In this Section, wo briefly recall how the IR diver- 
gences appear in the PT for an interacting boson system 
in the presence of a condensate beyond the Bogoliubov 
approximation, and then set up the RG treatment to 
heal these divergences. To this end, we will exploit a 
functional-integral approach which is especially suited to 
the purpose.^'' 



F[X, X*,A^] = In jDil>*Vil> g^I*''''*'^'^*'^-! . (2.2) 

In particular, the single-particle normal and anomalous 
Green's functions are given by 



Gn{xi,X2) = - {lp{xi)'tp*{X2))c = - 



6F 



6X*{xi)5X{x2) 

dF 



gMxuX2) ^ -m^^)^i^^))c = -j^^^^^sx*{x2) 

where the connected average {ip*tp)c = (V'*V') ~ (■'/'*)(V') 
is defined with respect to the weight function exp{5}. 

Bose condensation is conveniently introduced in terms 
of the Fourier transforms 



where k is a wave vector, ujg = 2TTsf3^^ (s integer) is a 
bosonic Matsubara frequency, and fl is the volume occu- 
pied by the system. We set 



(2.4) 



A. Bogoliubov approximation and the appesirance of 
IR divergences 

We consider a system of neutral bosons interacting 
through a short-range potential v, in the presence of aux- 
iliary external sources A, /i, and A that serve to generate 
the field, density, and current connected correlation func- 
tions. The associated action is given by: 

S = y"da;|v'*(x) [dr + fi{x)]tjj{x) - \{V - iA)il;{x)\^ 

-| mxt + i'{x)X*{x) + r{x)X{x)'j . (2.1) 

In this expression, x = (r, r) (with r ranging from to 
the inverse temperature (3) is a vector with d+1 dimen- 
sions, ip{x) is a bosonic field obeying periodic bound- 
ary conditions over the imaginary time r, and v is the 
two-body interaction. We shall always consider the zero 
temperature limit (3 ^ oo. Since we are concerned with 
the IR behavior, we can choose v local from the outset, 
with an ultraviolet (UV) cutoff implied in UV divergent 
integrals. We have also set h = 1, ks = ^, and m = 1/2. 

The sources X{x) and A(x) are let to vanish at the 
end of the calculation, while /i(.T) reduces to the ordi- 
nary chemical potential jj. In the following, we shall 
sometimes use the short-hand notation Ai, = {fi, A) with 
1/ = 0, . . . , d. 

The relevant correlation functions can be obtained 
from the free energy functional F, defined as: 



where now k = (a's,k) and a is arbitrary, and intro- 
duce Eq. (2.4) into the action (2.1). In this way, the 
action can be split into a (free) quadratic part Sq = 
Tlik V'*(^)^(^)S^^(^) with associated inverse propagator 



g-^{k) = i(jOs-'k^ + IJi 



(2.5) 



and into the remaining perturbation terms. ^ The param- 
eter a is then fixed by requiring that {tp{k = 0)) = 0, 
which corresponds to the absence of "tadpole" diagrams,^ 
i.e., to the vanishing of 1-particle irreducible diagrams 
with a single line entering: 



0. 



(2.6) 



The value of a thus becomes a function of the uniform 

external source A, and coincides with the (square root 
of the) condensate density in the physical limit A ^ 0. 
This requirement simplifies PT considerably. In particu- 
lar, it is convenient to introduce a matrix Green func- 
tion by defining ^i(fc) = i^{k) and V'2(fc) = V'*(~^)- 
In this way, Gnik) = Q22{-k) = - (i){k)'il)*{k)) and 



Qi2 = Qiii-k) = - {tl){k)i){-k)j satisfy the Dyson- 
Beliaev equation g{k)-'^ = G-'^{k) - T,{k) with 



r - f So{k) 
^"'y Goi-k) 



(2.7) 
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and the matrix elements of the self-energy E are distin- 
guished by the number of incoming and outgoing lines. 
The off-diagonal terms Si2(fc) = ^2i{k) characterize the 
broken-symmetry state as they vanish identically in the 
normal state. Given a simple matrix inversion 

leads to G- The Bogoliubov approximation consists in 
taking the diagrams with no loops for the self-energy: 
Sii(fc) = 2w|ap, T,i2{k) = Q?v. Within the same ap- 
proximation, the vanishing of tadpole diagrams is equiv- 
alent to the vanishing of the linear term in the action. 
This gives for A = 0: 



and the following form for the Green's functions; 



El 



(2.8) 



(2.9) 



2/xk^. Note that, within the Bogoliubov 



with El = \i'' 

approximation, the occurrence of the sound mode (i.e., 
i?k cx |k| for k^ <C 2y^) is strictly related to Si2(fc = 0) 
being non- vanishing. 

Inclusion of higher-order terms beyond the Bogoliubov 
approximation, however, spoils the above results due to 
the occurrence of IR divergences. For instance, let us 
consider the following one-loop diagram contributing to 
2^12: 



m 




2 2 r d 

a V Co 

J qHk 



(P{k + q) 



n\2 



(2.10) 



where A is an UV cutoff and q is a. {d + l)-dimensional 
vector. We have used the small fc-form of the Bogoliubov 
propagators (2.9) that behave behave like k~'^, where 
k = (ws/cojk) with = 2/i being the square of the 
sound velocity. Problems arise because A{k) diverges 
when < 3 as A; — > (in the zero-temperature limit) as 
follows: 



m 



k'^-^ for d < 3 

ln(fcVA2) ford = 3 



(2.11) 



while it remains finite for d > 3. The one-loop contribu- 
tion (2.10) thus overcomes the Bogoliubov contribution 
for small enough k. It is then natural to introduce in 
this context a generalized Ginzburg criterion for the mo- 
mentum variable, by comparing one-loop and zero-loop 
self-energies: 



Si°2^(fcG) 



1/2 



l,d-3 



d<3 



log{co/kG) d = 3. 



(2.12) 



Here, we have used w n and we set the UV cut- 
off A equal to Co as simple estimates. According to 



this criterion, IR divergences become dangerous when 
k < ko, where ka = {v^n) ^c^-d) for d < 3 and kc = 
Co exp{— (v^n)~^/^} for = 3 as given in Ref. 10. 

Notwithstanding the presence of these singularities, ex- 
plicit calculations have shown that IR divergences cancel 
out in most physical quantities.^ A notable example is 
the Hugenholtz-Pines identity^, namely. 



Sn(fc = 0)-Si2(fc = 0)=M 



(2.13) 



which implies cancellation of the IR divergences in the 

two self-energies. 

The identity (2.13) provides the simplest example of 
WI, whereby the (irreducible) diagrams of FT are con- 
strained by the underlying gauge symmetry. A system- 
atic study of the implications of WI was originally per- 
formed in Ref.'''. In the following Sections we will exploit 
extensively the whole set of WI to reduce the set of IR 
divergent diagrams of PT to only a few independent ones. 

As mentioned in the Introduction, not all IR diver- 
gences, however, disappear from the theory. In particu- 
lar, the longitudinal single-particle Green's function (to 
be discussed in detail in Section IV) has to diverge in the 
IR on physical grounds, implying that 



Ei2(A; = 0) = 



(2.14) 



for c? < 3. This was found as an exact result by NN. It is 
clear that the result (2.14) which depends on dimension- 
ality cannot be inferred from a WI like Eq. (2.13), which 
is instead independent from dimensionality. It is also 
clear that WI by themselves are not sufficient to specify 
completely the IR behavior of the system. To this pur- 
pose, specific dynamical properties have to be taken into 
account in addition to WI, as it will be discussed in Sec- 
tion IV where the RG equation for a single independent 
running coupling will be derived. 



B. Representation for the Field-Theory RG 
treatment of interacting bosons 

The RG treatment in the broken-symmetry phase is 
considerably simplified when replacing the Bose field V 

and ip* by their longitudinal (?/;;) and transverse {ipt) 
components to the direction of the broken symmetry: -"^^ 



ilp(x) = a + ij)i{x) + iijjt{x) 
ip*{x) = a + ^i{x) - iiptix) 



(2.15) 



where ipi and ipt are real functions, and the order param- 
eter has been taken real as well without loss of generality. 
These new fields are introduced to distinguish explicitly 
between transverse and longitudinal fluctuations, whose 
difference is not evidenced by the standard representa- 
tion resulting into Eq. (2.9). 
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In terms of Fourier transform of the new fields, the action (2.1) becomes for A = and fi{x) = ji: 



fe 

fcl+/C2 + fe3=0 



(2.16) 



fel+fe2 + fc3+fc4=0 



Quite generally, the order parameter is given by 

^U^) = (?A,(x)) = SF/6X,{x) (2.17) 

where i = {l,t), A = (A; + iXt)/2, and F is given by 
Eq. (2.2) (once expressed in terms of A; and A^). Con- 
sistently with our choice, we shall eventually assume the 
limit At = and A; 0, so that only tpi^ survives. 

The variable a can be eliminated by the tadpole con- 
dition. It can be useful to clarify this point in more de- 
tails here, since at higher orders elimination of a becomes 
more subtle. Given the action (2.16), one can regard the 
quadratic part as the free theory and treat the remainder 
as a perturbation, including the linear term in tpi (we set 
Ai = from the outset) . The diagrammatic perturbation 
theory is then well defined. Diagrams are conveniently 
summed up by grouping them according to the number 
of loops present in each of them. One can verify that this 
is equivalent to scale the action as S ^ S/a and calcu- 
late all quantities by an expansion in the parameter a. A 
naive way of proceeding would be to choose the value of 
a such that the linear term in the action vanishes, giv- 
ing V = fi. In this way, propagators are massless and 
all physical quantities can be calculated within the loop- 
expansi on w ith the parameter a fixed to its mean-field 
value \/Ji/v. Actually, the presence of the 3-legs vertices 
makes this approach cumbersome. As a matter of fact, 
all diagrams will have tadpole corrections. It is not diffi- 
cult to recognize that the effect of these corrections could 
be incorporated in a redefinition of a. The shift required 
would be exactly the correction to the average of tpi in- 
duced by the tadpole corrections: {tpi) = a + {ipi) ^ a. 
An alternative procedure is thus to fix a by the condi- 
tion that tadpole 1-particle irreducible diagrams vanish. 
This condition guarantees that (^() = a, since the van- 
ishing of 1-particle irreducible diagrams implies also the 
vanishing of all reducible ones contributing to (ipi)- The 
equation for a is then calculated perturbatively as a loop- 
expansion without tadpole insertions: 



ad) = 



where F(°)(a) = 2a(/i - va^) and F(")(a) is the 
sum of all one-leg n-loop diagrams (cfr. (Bl) and (B2) 
for the one- loop case). Substituting the loop expan- 
sion of a (a = a(°) + + . . .) and solving for 



a order by order, one obtains a^^^ = ii/v, 
-i^(i)(a("))/(aF(")(a)/9a)„^„(o), .... This procedure 
defines a much simpler diagrammatics, since no tadpole 
diagrams are present, their contribution appearing as a 
shift of a. The price to pay is that one has to calcu- 
late all diagrams for arbitrary a both in the 3-leg cou- 
plings and propagators. This implies that the propaga- 
tors would have a mass (gap), but since all quantities 
have to be expanded around their mean-field values one 
recovers massless propagators in all expressions.^^ This 
fixes the procedure in a precise way. 

In our case an important simplification occurs. As we 
are interested in the IR behavior of the exact propaga- 
tors, finite corrections are not important. Since the tad- 
pole diagrams are all finite in the infrared, we will not 
calculate the finite corrections to a accordingly. In the 
following, a can be then thought as fixed and to coincide 
with ilno. At the same time, the propagators are massless 
and no tadpole diagrams will be considered. In practice, 
this means that finite corrections to the condensate are 
not important down to = 1. At that point the con- 
densate vanishes due to IR divergences. We will see that 
precisely at d = 1 our approach breaks down since the the 
longitudinal and transverse fluctuations becomes equally 
singular. 

As far as the chemical potential entering the action 
(2.16) is concerned, it is possible to leave its value un- 
specified in the following calculation since the interact- 
ing boson system with a fixed chemical potential is well 
defined. 

PT simplifies by introducing the Legendre transform 
of F with respect to the sources A: 



r[iP,o,A,] = jdxX,{x}^Ux) - F[Xi[iPoi], A,] (2.19) 
where 

ST/Si^io{x) = Xi{x) (2.20) 



F^''\a)+F^'\a) + F^^\a) + ... = (2.18) and 



6r/SA^{x) = -6F/6A^{x) . 



(2.21) 



Similarly, the connected correlation functions and the 
vertex functions are respectively given by: 
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SXi^ {xi) ... 6Xi^ (x„)M^i (2/1) •• • <5^i>,„ {Vm) 

Qii...in;i'i...i^rn ('^l? • • • 5 '^n'l Ul: ■ • ■ : Urri) ; 



^(n+m)p 



'5V'iio(a;i) . . . Hir,o[xn)SA,,^{yi) . . . 6A^^{ym) 



(2.22) 



(2.23) 



Here, the derivatives with respect to give the con- 
nected and amputated (1-particle irreducible) correlation 
functions, for F and F, respectively The derivatives with 
respect to Ai, generate the current and density response 
functions (or their amputated counterparts). 

The single-particle Green's functions analogous to 
Eq. (2.9) can be obtained from the quadratic part of the 
action (2.16), which reads in terms of the independent 
fields: 



with 

M{k) 



Us 



-LOs 



jjL — va — 



(2.24) 



(2.25) 



With the aid of the condition (2.8), one then obtains: 



g{k) = -{M-')' 



1 1 / k2 



(2.26) 



which for small k and behaves like 



G^^{k) ~ -c2/(^2 + c2k2) 

Quik) ~ -C0,,/(W? + C2k2) 

QuOi) ~ -kV(u;2 + c2k2) 



In the present case the low-energy region is dominated 
by the Goldstone (sound) mode. The IR limit of the 
correlation functions is thus naturally defined by tak- 
ing frequency and wave vector to vanish while main- 
taining their ratio constant. For example, one can take 
lim~^-^Q Qij(jujs,j'k). In this way, the three propagators 
(2.27) have different IR behavior, with the strongest sin- 
gularity residing in the transverse propagator Qtt- This 
situation has to be contrasted with the standard tp- 
representation, where the Bogoliubov propagators (2.9) 
share instead the .same IR behavior. The choice of 
{ipi , ■^t) in the place of{ip,'ip*) will thus turn out to be cru- 
cial to select the interaction terms on the basis of their 
relevance. Furthermore, this choice will be important 
from a physical point of view since an anomalous dimen- 
sion will appear in the longitudinal component Gii only. 

To treat these singularities, it is important to classify 
the coupling terms in the action (2.16) according to their 
relevance in generating IR divergences. A standard way 



to perform this task is to define scaled fields and frequen- 
cies, so that the infrared dimensions of the propagators 
in (2.27) coincide with their engineering dimensions. In 
our case the scaling is defined as follows: 



V'((fco,k) = Co'ipi{coko,k) 
■0t(fco,k) = ■0t(cofco,k) 



(2.28) 



with the short-hand notation ko = icOg/co (from now on 
k = (fco,k)). Since the rescaled fields have dimensions 

Mk)] = [Mk]] + 1 = and [Mk)] = Mk)] = -1 
(we have defined the scale of the momentum as the unit, 
[k] = 1), this leads to the desired dimensions for the 
propagators of the rescaled fields: 

[Gttik)] = -2 , [Guik)] = -1 , [Giiik)] = . 

Since different diagrams contributing to the same quan- 
tity must have the same dimension, couplings of higher 
dimension must appear together with most divergent 
Green's functions. Thus the engineering dimensions of 
the couplings for the new fields can be used to classify 
their relevance. Equation (2.28) implies that, in real 
space, the fields are rescaled as 



— 1 /2 ~ 

t{jl{xo,x) = Co' 1pl{xo/Co,X) 



(2.29) 



which follows from a redefinition of the Fourier trans- 
form: 



VJJcJl Jo 



/3c„ 



dxo ld''xe-'''''ipi{x). (2.30) 



The dimension of the fields are [ipi{x)] = {d + l)/2 and 
[iptix)] — {d— l)/2. From the above considerations, the 
dimension of a generic coupling Wn,nt associated with the 



(2 27) monomial ■)/'"' (x)^/'"' (x) is given by 



, , d+1 d- 
Ki.nJ = -ni—^ rit — 



d+1. 



(2.31) 



For later convenience, we write explicitly the dimensions 
of the couplings up to four legs: 



[vt] 



e/2 
e/2 



M = -2+f 

[viit] = -1 + f 

M = f 



[vii] = 

[vit] = 1 , (2.32) 
[vtt] = 2 

[villi] = -4 + 6 

[viiit] = -3 + e 

[vim] = -2 + e (2.33) 

[vittt] = -l + e 

[vtttt] = e 



where vi = vio, vu = V20, vm = V12, etc., and e = 3 - d 
parameterizes the distance from the critical dimension 
dc = 3 at which logarithmic divergences appear. 

Couplings with more than four legs have not been con- 
sidered since they are irrelevant at d = dc, i.e., they 
have negative dimensions. The vertex functions have 



6 



the same dimensions of the relative couphng constants. 
This impHes, for instance, that the tadpole equation 
is non-divergent in the infrared, as anticipated, since 
[vi] = [Ti] = 2 - e/2 > 

To make PT calculations (as well as the RG treatment) 
more convenient, we generalize the original action (2.16) 



in the broken-symmetry phase as follows. Besides the 
terms already present in (2.16), wc introduce additional 
running couplings which are relevant or marginal at dc 
(while all irrelevant couplings were already present in the 
original action). The generalized action thus acquires the 
form:^^ 



{[vii + ziik^]i;i{-k)Mk) + 2ht + witLUs]M-k)Mk) + [vtt + uaujl + Ztt\<i^\i>t{-k)i>t{k)] 

k 



3! 



X! {vmtlJt{ki)tpt{k2)tlJt{k3) + Svitttpi{ki)ipt{k2)'tpt{k3) + vuiipi{ki)ilJi{k2)ipi{k3)} 

ki+k2+k3=0 



4! Y{vttttMki)Mk2)Mk3)Mk4) + Qvuui)i{ki)i>i{k2)Mk3)Mk4) + viui^lJii^^^ 



(2.34) 



The reason to use these new couplings in the place of 
the original ones is related to the fact that they appears 
as coefficients of the monomials of fields, and thus have 
a simple bare scaling behavior. It is also clear that there 
is a definite relation among the; new set of couplings and 
the original ones, that can be inferred by direct compar- 
ison of (2.16) and (2.34). Since the new set of couplings 
is larger, it is also clear that they will not be indepen- 
dent when related to the original parameters. In Section 
III A these connections will be discussed and related to 
the gauge invariancc of the original theory. 

As far as the marginal couplings (wu, utt, zu) arc con- 
cerned, we recall that in the original action (2.16) they 
had, respectively, the values (2,0,2) regardless of the orig- 
inal parameters (w, /x, a). We follow here a standard pro- 
cedure of the field-theoretical RG treatment and let the 
coefficients of the marginal terms proportional to u), u"^, 
and assume arbitrary values in order to renormalize 
the theory. In principle, an alternative procedure would 
be to introducing two wave-function-like renormalization 
parameters and one dynamical renormalization param- 
eter which, in turn, allows one to keep the matrix of 
the propagators unchanged. It will turn out, however, 
that the renormalization of these parameters, although 
nontrivial, will not affect the scaling form of the matrix 
of the renormalized propagators, since these parameters 
will appear there only in some specific combinations that 
are independent of the RG fiow. It will, in fact, turn out 
that the anomalous behavior of the longitudinal propaga- 
tor depends only on the renormalization of the coupling 
vii; at the same time, the dynamical index will stay finite 
at 2; = 1, consistently with the Bogoliubov sound mode. 

Quite generally, the matrix entering the quadratic part 
of the action (2.34) is defined as in (2.24) and is given 
by: 



M = - 



vii 



ziik"^, 



Vlt + WltLOs 



Vtl - WltUJs, Vtt + UttO^s + -^ttk 



(2.35) 



in the place of Eq. (2.25). The corresponding single- 
particle propagator becomes: 



g{k) = - 



Vtt + Utiles + ^tt^^ 
-Vlt - WltUJs , 

Qu{k) Gitik) 
Qti{k) Gttik) 



-Vtl + WitUJa 
Vll + Zll^S? 



1 



D{k) 
(2.36) 



with 



D{k) = vaVtt - vft + {uttVn + wfjw^ 
+ {viiztt + vttzii + ziizttk'^)k^ 



(2.37) 



in the place of Eq. (2.26). For later convenience, the com- 
ponents of the propagator (2.36) are depicted schemati- 
cally as follows: 



Giiik) 



Qit{k)= , Gttik) =■ 



(2.38) 



The remaining terms of the action beyond the quadratic 
part are considered as perturbation. The interaction 
terms retained in the RG treatment are: 



vitt = ■ 



Vtttt 



-X 



(2.39) 



PT treatment then proceeds now along the usual lines. 

Although the generalized action (2.34) contains four 
marginal couplings [vu, wu, Utt, ztt) and seven relevant 
couplings (w;, vt, Vlt, Vtt, vm-, viu, vtut), it will be shown 
in the next Section via the use of WI that these couplings 
are actually not independent from each other. 
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III. GAUGE INVARIANCE AND WARD 
IDENTITIES 

The original action (2.1) satisfies local-gauge invari- 
ance, in the sense that it is unaffected by the following 
transformation: 

X{x) A(a;)e-*'^(^) X{x)* X{x)*e'^^''^ 
nix) n{x) - idrx{x) A{x) A{x) - Vx{x) ■ 

(3.1) 

This invariance provides definite connections (in the form 
of WI) among the vertex functions (2.23).^^^ These iden- 
tities have been exploited by Gavoret and Nozieres to 
establish connections between vertex functions and ther- 
modynamic derivatives for vanishing wave vector and fre- 
quency. In this Section, we shall make use of the WI in 
two additional ways. On the one hand, we shall con- 
straint the various couplings of the generalized action 
(2.34) by showing that they are not independent from 
each other when the generalized action is required to be 
gauge invariant. On the other hand, we will exploit the 
WI to connect the IR singular parts of the vertex func- 
tions in PT and to implement their cancellation to all 
orders in the response functions. 

To be consistent with our notation, we first give a short 
derivation of the set of WI relevant to our purposes. In 
terms of the longitudinal and transverse components {tpi , 
ipt), the original action (2.1) has the following invariance: 



S[Riji)j, RijXj, A, - d^x] = 5[V'i, Ai, A,] (3.2) 

where i and j take the values {I, t) and a summation over 
repeated indices is implied. In Eq. (3.2) we have intro- 
duced the notation = {idr,V) and R is the rotation 
matrix 



R 



cosx(a;) -sinx(a;) 
sinx(a;) cosx(a;) 



(3.3) 



x{x) being an arbitrary real function. The above invari- 
ance leads, in turn, to the following invariance of the 
functional F: 



T[Av - dt,x{x), Riji/Jjo] = r[A^, Vi. 



(3.4) 



The desired connections among the vertex functions are 
obtained at this point by differentiating Eq. (3.4) with 
respect to x- One obtains: 

Ti{x)aijiljjo{x) + d,T,,{x) = (3.5) 

where a is the infinitesimal generator of rotations 



-1 

1 



(3.6) 



and X has been set to zero eventually. Successive differ- 
entiations with respect to ipio give in addition: 



^im.{xi,X2)(Tij'll;jo{xi) +Ti{xi)aimS{xi - X2) + d^^Tm;,,{X2;Xi) = 0, 



(3.7) 



and 



^imn{xi, X2, X3)aij1pjo{xi) + Tim{xi, X2)crinS{xi - X3) + Tin{xi, X3)aimS{xi - X2) 
+d^^^mn;i^{x2, X3; Xi) = , 



■•{xi,X2,X3,X4)crijll)jo{xi) +Timn{xi,X2,X3)crir6{xi - X4) + Timr{xi, X2, X4)ain5{xi - X3) 



+Tinr{xi,X3,X4)(JimS{xi - X2) + {X2,X3,Xi\Xi) =0. 



(3.6 



(3.9) 



The WI (3.5)-(3.9) are sufficient for the RG treatment 
to be discussed in Section IV. 

For homogeneous external sources, in terms of Fourier 
transforms (for instance, Ti{x) = ^ e*'^^Fj((/) with 
the short-hand notation qx = u>sT + k • r) the WI (3.5)- 
(3.9) simplify considerably. In particular, the sets of WI 
which will be explicitly used in the following for vanish- 
ing transverse source (i.e., Vto = 0) can be grouped as 
follows: 

rtV/o = o, (3.10) 



rti{k)^io + Tt - iKTi,,{-k) = (3.11) 
Ttt{k)i^io -Ti- ik,TfA-k) = (3.12) 



from (3.5) and (3.7), which encompass the Hugenholtz- 
Pines^ identity for k = 0; 

^tti{ki, k2)'ipio + Ttt{-k2) - Tii{ki + ^2) 
-i{ki)^Tu,u{k2, -fci - fcs) = , (3.13) 

Ttttiki, k2)^io - TH{-k2) - Titiki + k2) 
-i{ki)uTtt.4k2, -fci - fcs) = (3.14) 

which follow from (3.8) and are the standard WI associ- 
ated with the continuity equation modified by the pres- 
ence of the 3-legs vertices; and from (3.9) the four legs 
ones 

^tttt{ki,k2, k3)ipio - Titt{-k2 - ks, ^2) 
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-^ittiki + h, k2) - Vitt{ki + k2, fca) 
-i{ki)^Tttt;i^{k2,k3, -ki - - fcs) = . 

(3.15) 

Here, k„ = {iug,^) with v = (0, . . . ,d) and the conser- 
vation of momentum has been used to eliminate the de- 
pendence on the last momentum of each vertex function, 
i.e., Ti{k) = Tiik = 0) = Ti, Tij{ki,k2) = Tijiki), and 
so on. It is thus understood that the vertex functions are 
free from the overall momentum conserving delta func- 
tion. We include the i in the definition of fco only in the 
scalar product entering the WI to simplify the notation; 
everywhere else ko = LUg. The rescaling (2.28) has not 
been considered explicitly in the present Section, but it 
can be readily introduced whenever necessary. 

We recah that the WI (3.10)-(3.15) have been ob- 
tained using the invariance properties of the original ac- 
tion (2.16). For arbitrary values of the couplings, the 
generalized action (2.34) will not be invariant under the 
transformation (3.2). Thus the requirement of fulfilling 
(3.2) will provide constraints on the couplings appearing 
in (2.34). It is also clear that this requirement by itself is 
not sufficient to reduce the number of independent cou- 
plings to the original ones. It is, in fact, possible to add 
to the original action (2.1) terms of the type w„|V-'P" still 
preserving gauge invariance. To identify (2.16) exactly 
with (2.34) one should specify, for instance, that no 5-leg 
couplings arc present. With this additional condition, it 
can be shown that the independent couplings left are only 
four (for instance, vi, Vu, ipio, and 4'to) corresponding to 
/i, V, Xi, and At of the original action (2.1). 

In Appendix A, we discuss in more details the formal 
definition of gauge invariance in terms of the couplings 
introduced in this Section. There we show how to en- 
force the invariance (3.1) directly on the action. In an 
equivalent way here we proceed by assuming that the 
new action (2.34) shares the same invariance properties 
of (2.1). Then Eq. (3.4) holds and the above WI follow. 
We thus enforce the vertex functions associated with the 
new action to satisfy these WI and obtain relationships 
among the bare couplings of the action (2.16). 

To be more precise, we can calculate all vertex func- 
tions by a loop expansion, since the WI are preserved 
order by order. We can use this method to: (i) Con- 
straint the bare couplings using the lowest order 0-loops 
expressions; (ii) Connect the divergent contributions to 
different renormalized running couplings at any loop and 
thus find relations among themselves. 



A. Global Ward identities: k = 

We are now in a position to discuss the vanishing of 
the relevant bare couplings vi,Vt, vn, Vtt, and Vut present 
in (2.34) at any dimension. 

At the lowest order in the loop expansion, Eqs. (3.10)- 
(3.12) and (3.14) can be evaluated at vanishing k and the 



vertex functions expressed in terms of the bare couplings. 
This procedure provides definite relationships among the 
relevant bare couplings. Specifically, we obtain vt = 
from Eq. (3.10) for non vanishing condensate density tpio', 
Vit4>io + Vt = from Eq. (3.11); Vuipio - Vi = from 
Eq. (3.12); and Vttti'lo - ^vu = from Eq. (3.14). This 
gives: 

vt = 0, vit=0, Vttt = 0, (3.16) 

and 

vi = tpioVtt (3.17) 

for vanishing i/'to- To determine vi, we consider in ad- 
dition the "equation of state" (2.20) (with i = I) for an 
isolated system with Aj = 0. This is the counterpart of 
the condensate equation (2.18) in the new variables and 
fixes the bare coupling vi = 0. Actually, this is not true 
at higher orders, since the shift of the condensate will in- 
duce finite corrections to vi that must be fixed in order to 
let tadpole diagrams vanish. Since we are not interested 
in finite corrections, we can set the bare and renormal- 
ized couplings = 0, and thus Vu ~ The vanishing of 
Vtt implements the Hugenholtz-Pines identity, leading to 
gapless propagators. 

We pass now to discuss the connections among the 
marginal couplings. Let's consider the hare marginal cou- 
plings first. The identities (3.13) and (3.15) can be in- 
terpreted at the lowest-order of the loop expansion as 
follows: 

vitti'io -vu=0 (3.18) 

VtttttlJio - ivitt = (3.19) 

where all k have been taken to vanish. The marginal 
bare couplings thus eventually reduce to four indepen- 
dent ones vii, Wit, utt, and ztt- 

It is clear that the relations (3.16), (3.17), (3.18), and 
(3.19) must be enforced on the bare running couplings to 
guarantee the gauge invariance of the theory. A field the- 
oretic RG treatment is now possible only if we neglect all 
irrelevant running couplings. Although this spoils gauge 
invariance, it does not matter for what concerns the lead- 
ing order-singularities. We can verify that gauge invari- 
ance is lost if we set all irrelevant running couplings to 
zero, since in this case the WI (at the lowest-order in 
loops) are no longer compatible with finite values of the 
marginal running couplings. In other words, the use of 
a perturbation theory for the F's with only the marginal 
running couplings will violate the WI (as shown explicitly 
in Appendix B), but the leading diverging contribution 
to each T function is fully retained by this approxima- 
tion. This implies that the finite part of the WI can 
be violated, but the divergent contributions coming from 
different T functions will correctly simplify in the WI. 
Thus the WI hold for the renormalization of each vertex 
function, since these are introduced exactly to subtract 
out the divergent parts. For this reason the WI (3.18)- 
(3.19) will hold for the renormalized running couplings 
too. 
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B. Local Ward identities: small k 

There remains to consider the connections between the 
couplings [wu, uit, zu), associated with the wave vector 
and frequency dependence in the action (2.34), with the 
thermodynamic derivatives in the limit of vanishing k 
and uj. To this end, we consider the vertex functions F;;, 
Tit, and Ttt associated with the couplings of interest. The 
long- wavelength limit of these vertex functions has been 
obtained long ago.''"^ Here, we will rederive them in a 
more compact form taking into account the presence of 
IR divergences. Irnporta.nt information on the scaling of 
the running coupling will follow from these connections. 

We begin by establishing their symmetry properties. 
From previous considerations Tit{k = 0) = and Tii{k = 
0) = 0. Using the symmetry properties of the action 
(2.34) under space and time reversal and under exchange 
of the two components I and t of the bosonic field, one can 
show that F;t {k) is odd under the replacement ujs —ujs 
while Ttt is even. We then take the zero-temperature 
limit (whereby Us becomes a continuous variable, de- 
noted by LO from now on) and consider the WI (3.11) 
and (3.12). 

It is convenient to cast the vertex functions in a form 
that makes their symmetry properties manifest. For in- 
stance, from r;t(a;, k) = — r;t(— a;,k) we write 



Tit{oj,-k)=ujVit{ujX) 



(3.20) 



where VuioJ, k^) is an even function of Wg- Similarly, we 
can write 



and 



Ti.o{k)=Vi.,o{ujX 



Tl;n{k) = U}iknVl;v{l^,'k^) 



(3.21) 



(3.22) 



(n = 1, . . . ,d) where Vifi{k) and Vi-v{k) are even func- 
tions of u) and depends on k only through k^ (note that 
the suffix V is here a reminder that we are considering a 
vector function). In this way, the WI (3.11) becomes: 



i^ioVitik) + Vi-oik) + k^ViAk) = 



(3.23) 



where Eq. (3.10) has been taken into account. 

To the purpose of understanding the behavior of the 

singularities when we let or k vanish indepentently, 
we discuss the form of the divergences. At the critical 
dimension {dc = 3), a generic vertex fmiction F(A;) of a 
single external variable in the limit k ^ can be written 
formally as 



T{k) = (fc2)^/2 
+ C 



(3.24) 



where D stands for the IR dimension of F(A;) and Am, 
B, and C are constants. [In general, instead of the term 



j-^2-j_D/2 Qj^^ ^jgg ^^2-^D/2^ ^^^2-^0/2^ ^j. ^^^^ linear com- 
bination of the two.] In this expression, we are not con- 
sidering the terms arising from the internally diverging 
diagrams, since the divergent contributions due to these 
diagrams are eliminated by the renormalization at lower 
order of the other T functions. In this scheme, the di- 
vergences left are all due to the primitively divergent di- 
agrams of each vertex function. These divergences must 
cancel out in the WI, so exact relations result among the 
primitively divergent contributions to each F function. 
The form (3.24) implies also that the fc — > limit of the 
divergent log series does not depend on which limit is 
taken first, w ^ or k ^ 0. This is a consequence of the 
fact that all singular dependence comes through the sin- 
gle argument of the log terms. This fact can be verified 
order by order in loops [see Appendix B] . Of course, the 
order in which the limits are taken affects the final value 
of the F functions that have overall u) and k terms. The 
F functions that have D = will be totally symmetric, 
and the information gained for one of the two variables 
set to zero applies also for fc ^ 0. Knowing the degree of 
divergence of the vertex functions and the explicit com- 
binations of UJ and k appearing in each expression, one 
can establish relations between the vertex functions in 
presence of IR divergences. 

In order to proceed further we need the IR dimen- 
sion (i.e., the engineering dimension after the rescaling 
(2.28)) of the composite vertex parts. One should con- 
sider, however, that IR dimension is correctly defined for 
monomials of fields. For instance, we can readily obtain 
that the dimension of r;.;2 is 2 — e/2, where by l'^ we mean 
the field tjjf = — (tpf{x))). But the dimension of 

a physical composite fields will mix different monomials 
and thus different IR behavior: F^^o = 1";:^ +r;./2. As 
we are concerned with divergences, we identify the IR 
dimension of a composite field with the most divergent 
dimension (in the above example with [F;.t2] = —e/2 we 
define [F,;o] = -e/2). 

We report the IR dimension of other relevant F func- 
tions for further reference: 





= 2 




r;00 




= 1 




r;iO 


ft;/ 


= 2 




F 



— e 

2-e 

4-e. 



(3.25) 



These identities are obtained by performing the Fourier 
transform of the Green's function, taking into account 
that square of fields depend on a single variable. 

We pass now to discuss the WI (3.23). Since the 
functions entering (3.23) have vanishing IR dimensions 
[cf. Eqs. (3.25) and (2.33) taking also into accont that 
the functions V differ from the respective F by factors to 
or k], the last term on the left-hand side of (3.24) van- 
ishes as k^ ln'"(k^ -|- w^/c^ in the IR limit. In this way 
we obtain: 



■^io lim Vit{k) = - lim Vi-Q{k) 



(3.26) 
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(which is related to wu) where 



dT 



limP;;o(fc) = — 

fe->0 pil OlploOIJL 



(3.27) 



is a thermodynamic derivative by definition. We can also 
obtain the same results by setting k = and letting 
w ^ 0. The final result is valid for k 0, since Vi-Q{k) 
is a function of + w^/c^ only for fc ^ 0. 

We consider next the WI (3.12). Expressing again 
the vertex functions V in terms of the V functions, 
r^.Q = wPt;0 and Tt-i = ikiPt,v as in Eqs. (3.20)-(3.22), 
we write: 



from which we obtain for 



lim 



rtt(o,k) 

k2 



-^limPt;,(0,k2) 



(related to Zu), and for k = 



lim 

OJ^O 



r«(w,o) 



-3- lim Pt;o(w,0) 



(3.28) 



(3.29) 



(3.30) 



(related to Uu)- Note that the function Vt-.v has dimen- 
sion 1 — e/2. This suggests that the primitively divergent 
contribution to this function vanishes at dc- By inspec- 
tion of the PT one can verify that this is actually the 
case since the current insertions reduce the strength of 
the divergence. This fact is important since it implies the 
absence of singularities on the left-hand side of Eq. (3.29) 
(and thus the fact that will not scale). This docs not 
imply that Vt-v vanishes in the limit of fc ^ 0, since the 
non-singular contribution can be finite. 

The function Vuv^k) in the limit fc ^ cannot be re- 
lated directly to thermodynamic derivatives, since these 
involve derivatives with respect to the transverse field 
V'to- In order to identify the two limits (3.29) and (3.30), 
we use the following additional WI (obtained by differ- 
entiating Eq. (3.5) with respect to and then taking 
the Fourier transform): 



(3.31) 



which in terms of the V functions reads (with a factor uj 
divided out) 



i^ioVffiik) - V,oo{k) - k2p.„o(fc) = 



(3.32) 



and 



ikni^ioVt-vik) - iJ^knV-ovik) 



k2 



h h 
~k2~ 



(3.33) 



(to, n = 1, . . . ,d), where a sum of repeated indices is im- 
plied and the quantity within brackets is T-mn- Equation 
(3.32) for k = gives 



tpio lim Vffiik) = lim T-ooik) 



1 d^r 



(3.34) 



which is again a thermodynamic derivative and relates 
Utt to Fjoo- Multiplying Eq. (3.33) by — and sum- 
ming over n we obtain further 



i^ioVt-v{k) - w^^O;0^(fc) -h 0;(fc) = , 
which for a; = and k ^ yields 



(3.35) 



'4)10 lim Vt;v(k) = - lim Mk) = . (3.36) 

k^O k->0 m 

Here, m{— 1/2) is the mass and the last identity fol- 
lows from the Galilean invariance of the action.^ Equa- 
tion (3.36) relates Zu to n/m. Finally, in the limit fc — > 
the limit oiTu{k) is {PVL)-^ d'^V / dip^^ by definition. 

We have thus verified that, even in presence of IR 
divergences, the asymptotic limit of the quadratic ver- 
tex functions and of its derivatives arc related to ther- 
modynamic derivatives. Summarizing, in the IR limit 
one expects that Tit{k) = —ri-o{0)u/ijjio and Ttt{k) = 

Using the scaling of the RG obtained in the next Sec- 
tion, we obtain further that the quantities 

Til (pk) ~ vii (p) , Tit (pk) ~ wp Wit (p) , 

Tttipk) ~ y u„{p) + k^ zttip)] (3.37) 

when the RG scale p (defined in the next section) goes 
to zero. 

These connections between the k ^ limits of the ver- 
tex functions and the asymptotic limits of the RG cou- 
plings allow us to identify the asymptotic p ^ limit of 
the ratio 



Viiip) Zttip) 



vii{p)utt{p) + witip)^ 



(3.38) 



with the macroscopic sound velocity Cg- To this end, we 
write 



lim c^ip) = 



g2p j^n 



n/m 
(dn/dp)} 



1 d^T 



i^fo \di^iodp 



(3.39) 



which equals the square of the macroscopic sound veloc- 
ity by standard thermodynamic arguments. This last 
identification is obtained by expressing derivatives at 
constant tpio in terms of physical derivatives at constant 
external sources A. The ratio (3.38) enters the determi- 
nant (2.37), where it identifies the (square of the) micro- 
scopic sound velocity. The constant term {vuVtt — vf^) 
vanishes identically as explained above. 

Finally, we note in this context that the ratio 
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'lim^ 

' p^O Wlt{p) 



Wo 



1 d'r 



dtpio 



(3.40) 



is proportional to (the inverse of) the "condensate" com- 
pressibiUty taken at constant external sources, which 
makes it a truly physical quantity. We thus expect on 
physical grounds that the ordinary compressibility dn/d^ 
and the "condensate" compressibility duo/d^ (as well as 
the density n) are all finite, since the physical system we 
are considering is in a stable phase. As a consequence, 
singular PT contributions will have to cancel out exactly 
in these particular combinations of running couplings, so 
that the ratios (3.38) and (3.40) as well as Ztt{p) (which 
would give an anomalous dimension to the transverse 
field) will not change under the RG flow.-'^'^ In Section 
IV we will present a diagrammatic proof of this fact for 
the two couplings Uu and wu, while the argument given 
after Eq. (3.30) is valid to all orders in PT. 



IV. FIELD-THEORETIC RG TREATMENT 

In this Section, we set up and solve the RG equations 
for the flow of the coupling constants discussed in the 
previous Sections. We will explicitly calculate the one- 
loop PT contributions to the relevant vertex functions, 
and show that higher orders of PT do not modify the 
IR behavior found at the one-loop level. In other words, 
we will be able to assess the exact IR behavior of the 
interacting Bose gas at zero temperature. This result 
contrasts with the standard (e.g., (p'^) theory of critical 
phenomena, whereby critical exponents can be calculated 
only approximately. The physical origin of this simple 
behavior in our case stems on the fact that the same 
theory is free from IR divergences when formulated in 
the variables phase and amplitude. In our variables the 
theory remains interacting at the fixed point to give a di- 
verging longitudinal susceptibility. Nevertheless, we will 
show that the effective interaction entering the transverse 
vertex part Ttt vanishes exactly. This explains the dif- 
ference with the (j)'^ theory in the final result for the Gtt 
correlation functions and recovers in an effective way the 
asymptotic freedom for the phase-amplitude action. 

As usual, we shall discard all irrelevant running cou- 
plings of (2.34) and write equations for the four marginal 
couplings {vii,wit,Utt, Ztt)- Before considering our spe- 
cific problem, we summarize the RG procedure employed 
shortly.23.14 



A. RG procedure 

Given the set of bare nmning couplings {g°} of the the- 
ory (regularized in the UV with a cutoff A), we impose 
the following normalization conditions: 



rj(A;)|fe2=^2 =gi. 



(4.1) 



The renormalized vertices are defined in the following 
way: 



^R{kj,9ii9°, K, A),k)= T{kj , gl A) . 



(4.2) 



Differentiating Eq. (4.2) with respect to k at constant g° 

we have: 



d_ 

dn 



d_ 

dgi 



Tii{ki,g,K) = Q. (4.3) 



If the dimension of the coupling gi is ai, we can introduce 
the dimensionless couplings Ui and u°: gi = n'^^Ui and 
g° = k"*u". The RG equation takes then its standard 
form 



d_ 

dn 



where Pi{u) = K{dui/dK)go. 
the following equation holds: 



TR{kj,Ui,K)=0, (4.4) 
One can then verify that 



Tnipki, u, k) = p^Tu^ki, u{p), n) , 



(4.5) 



where p dui{p) / dp = Pi{u) and a is the bare dimension of 
the vertex function. To obtain the IR behavior of a ver- 
tex function, it suffices to scale it with its bare dimension 
a and to substitute the new value of the running coupling 
u{p). In particular, when the normalization condition is 
^R{k, K)|fe=„ = 9i = K"'Ui(l) we obtain: 



TR{pki, u, k) = K"*p"'Ui{p) . 



(4.6) 



In this way, if near the fixed point the running coupling 
is flowing to zero with an exponent yi [ui{p) ~ p^'], the 
vertex part will have the following IR behavior: 



(4.7) 



It is thus simple to convert the p-behavior of the run- 
ning couplings into the fc-behavior of the relative vertex 
functions. The simple relation between the renormalized 
couplings and the vertex functions has been used to iden- 
tify the A; — > limit of Tij{k) given before (3.37). 



B. Additional constraints on the running couplings 

We now proceed to write four RG equations for the left 
four couplings: vu, wu, uu, and zn- All other couplings 
can be readily obtained from them. Indeed, the conden- 
sate ipio entering the WI is the exact condensate density 
which is not scaling. The other marginal couplings vm 
and Vtttt have been eliminated in favour of vu through 
Eqs. (3.18) and (3.19). 

It turns out that also the four left nmning couplings 
left are not renormalizing independently. As it was 
argued at the end of Section IIIB, we expect that 
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the combinations of couplings zti, wu/vu, and c^{p) = 
vuZtt/ {viiUtt + wfj arc invariant of the RG flow. Here, 
we show how these properties for wit/vu and c can be 
proved by analyzing the PT. These identities are shown 
to hold at one-loop level in Appendix C. The coupling 
Ztt has to remain c:onstant for the argument given after 
Eq. (3.30). The absence of divergences in this case fol- 
lows directly from the WI and power counting. We now 
proceed to prove the other two identities by analysis of 
the PT to all orders. 

The coupling wu can be identified with vu via the WI 
(3.13) and (3.14), which relate wu to r(/:o and vu to Tut, 
respectively. By inspection of the leading singular terms 
to all orders in perturbation theory, Ttt-o and Fjtt are then 
found to be proportional to each other. As a matter of 
fact, the only way in which we can make an insertion 
of the longitudinal (Z) leg is to use the bare interaction 
term vut, since it connects the external longitudinal leg 
to two internal transverse legs. It is then clear that this 
corresponds to performing a density insertion with the 
most relevant term (since \tp\'^ = tp^ + tpf + 2tpiotpt + "fpfo) ^ 
apart from the value of the bare external insertion that 
is different from the bare vm- (Sec Appendix A for a 
discussion on the density fields and the associated cou- 
plings.) Diagrammatically this reads: 



^tt;0 = 




ttl 



(4.8) 

We can write the chain of identities in the following way 
(by ~ we mean the leading-order behavior as fc ^ 0); 

i^Mko) ^-'^ r,„(o,fco) ^""-^ ^TttAo,ko) 



7tt;0 



WI vitt r/f(fco) 



7«;0 ko 



(4.9) 



where by •ytt-fi we indicate the perturbative insertion cor- 
responding to Ttf.o that in our case is wn (see Appendix 
A). In Eq. (4.9) we have set the spatial part to zero ex- 
ploiting the hypothesis of singular dependence through 
only Aiq -h k^. We thus conclude that, at the leading or- 
der in the divergent contributions: 



Tuiko) 



Wit ko 



(4.10) 



In this way, the logarithmic divergences entering in the 
RG equation for vu must also appear in the RG equation 
for Wit- If the exact (to all orders) RG equation for vu 
has the following form 



dvii 

P-^ = f{vu,wit,utt) . 



(4.11) 



the resulting equation for wu will be after a transient 
region: 



dwit , .Wit 
p—r- = j{vii,wit,utt) — . 
dp vii 

Dividing Eq. (4.11) by Eq. (4.12) we thus get: 

dwit Wit 



dvii vu 



impling 



Mp) 

wit{p) 



constant = Ci . 



(4.12) 



(4.13) 



(4.14) 



In a similar way, the invariance of c(p) follows from the 
exact connection between the singular parts of r;oo and 
Til, associated respectively with utt and vu: 



\ I 

A.lUM ^ Vltt 



r;oo = 



V7«;0 



(4.15) 



This procedure leads to the following asymptotic equa- 
tion for Utt 



dutt 

P-i— = -j{vii,wit,utt)—. 
dp v. 



(4.16) 



Using Eq. (4.14) we obtain dutt/dvu = — ^ and 

Utt{p) = -C^^vu{p) + C2, (4.17) 

so that uii + wfi/vii = C2 and c(p) is invariant (cf. 
Eq. (3.38)). The values of the two constants Ci and 
C2 can be extracted from the p — > limit of the run- 
ning couplings, specifically from Eqs. (3.39) and (3.40). 
The flow of the four running couplings will satisfy thus 
these three identities. We are left eventually with only 
one independent running coupling, for instance, Vu. 

The constraints found above for vu, wu, Ztt, and Utt fix 
completely the IR behavior of two out of three Green's 
functions, say Git{k) and Qtt{k). From Eq. (2.36) one can 
in fact say that the expressions for Git{k) and Gtt{k) at 
vanishing vn and vn involve only the above combinations 
of running couplings, which do not scale and are related 
to physical quantities. The resulting exact IR behavior 
is thus given by: 



att(k,a;) 



c mno 



1 



mc^ duo 



2n dp. ufl -\- c^k^ ' 



(4.18) 
(4.19) 
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This result is valid at all orders of perturbation theory 
and does not depend on the sealing of vu . The Goldstone 
mode is accordingly stabilized only by symmetry require- 
ments, irrespective of the behavior of vu. The vanishing 
of Vu (~ 1112) is no longer changing the IR behavior of Gtt 
and Git- As discussed in the Introduction, the IR behav- 
ior of Qii is expected to depend on dimensionality, and 
thus cannot be fixed by symmetry requirements only. 



C. One-loop equation for vu 

To proceed further we need to perform a specific calcu- 
lation of the RG equation for one coupling, for instance 
Vu. At the one- loop level, only one diagram contributes 
to Tu (see Appendix B). (Derivation of the other one- 
loop equations is not required. For comparison, we give 
the full set of one-loop equations in Appendix C for the 
marginal running couplings entering the propagators.) 
With dimensional regularization and after the e expan- 
sion has been carried out, Tu reads: 



T^u{k) = Vu - 



VuK 



2 ^2 , 



(4.20) 



where = k^ + fc^ and the factor = 2/{A'KY''^V{d/2) 
has been absorbed in the definition of the running cou- 
pling. We made use of Eq. (3.18) to connect vut to vu. To 
set up the RG equations, we need to define dimensionless 
running couplings (see Section IV A). Since all quadratic 
running couplings are already dimensionless (at any di- 
mension), only the condensate density V/o needs to be 
scaled by a factor k^/"^: 



e/2 



(4.21) 



Note that the powers of Cq have already been eliminated 

3/2 

in tjjio which is redefined like Co V-'Zo- By minimal sub- 
traction of the 1/e pole one obtains the following RG 
equation for vu{p)- 



dvu 1 



2 ' 



(4.22) 



where i^Hoip) = ipio{^)p'^^'^ exactly while ztt is invariant. 
We stress that the scaling of ipio is a trivial effect of its di- 
mension and the condensate density is not renormalizing. 
An alternative procedure would be to fix the condensate 
density and use scale dependent WI in order to connect 
the running couplings Vu, vm, and Vutt- The solution of 
Eq. (4.22) is straightforward and depends on dimension- 
ality. For d = 3, ipio is constant and 



vu{p) 
vu{l) 



1 -1 



Inp 



(4.23) 



For d < 3 we obtain instead: 



Mp) 

vu{l) 



1 + 



vu{l) 



2V'io(l)2z«(l)2 



(4.24) 



We recall that the RG equations are valid only after 
a transient region, so that the boundary conditions for 
Eq. (4.22) should involve the running couplings at some 
intermediate value < p < 1. For simplicity, we write 
the solution of this equation in terms of the coupling at 
p = 1. Since the results do not depend on the boundary 
conditions (provided the sign is not changed), we assume 
that the p = 1 is the point where the RG equations start 
being accurate. 

In both cases (4.23) and (4.24) the coupHng vu is flow- 
ing to zero, implying the vanishing of Si2(fc ^ 0) (here 
S12 = [Tu - Ttt + iiTit + Tti)]/^). At the same time, the 
constraint on the other couplings implies that Wu (p) van- 
ishes like vu{p), which from Eq. (3.38) gives that uuip) 
flows to Ztt{p)c'^{p)/Cg. Since both c{p) and Ztt{p) are 
constants of motion of the RG flow, Utt{p) — > -2«(1)- We 
can thus find the asymptotic expression for the longitu- 
dinal Green's functions as follows: 



Gu{k) 



1 utt kl + ztt^^ 

VllUtt kl + (c2/c2)k2 



Vu 



According to the scaling (4.7) of the RG, this implies that 



Gu{pk) 



vu{p) ' 



(4.26) 



The above result, together with Eqs. (4.18) and 

(4.19), gives the complete IR behavior for all one-particle 
Green's functions. By using the results of Section III on 
the finite value of the small /c-limit of the vertex func- 
tions, it is possible to rewrite the result (4.26) in terms 
of thermodynamical quantities: 



Gu{k,oj) 



8nnn 



dp 



CO 



r-.-r-i^k/K) 



(rf = 3) 
(1< rf< 3) 



(4.27) 



where A is an ultraviolet cutoff of the order of the sound 
velocity (corresponding to the region where the RG starts 
to be accurate) and Co is the sound velocity at the scale 
A. After that point, c{p) is, to a good accouracy, an 
invariant of the RG flow, so that we can substitute the 
fixed point value Co = c. 

We note that for d = 1 (e = 2) the transverse and 
longitudinal divergences become identical while the con- 
densate vanishes. It is thus expected that transverse and 
longitudinal fluctuations become symmetric since it is no 
longer possible even to define longitudinal and transverse 
directions. The main assumptions of our approach will 
thus break down, indicating the existence of a lower crit- 
ical dimension. Nevertheless, since the results arc valid 
exactly for 1 < d < 3, as said already, the predicted IR 
behavior approaches the correct result for d — > 1.^ 
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The usual Green's functions can be readily obtained 
from Eqs. (4.18), (4.19) and (4.27): 

^ c^muo 1 mc^ duo ito 

+guiK'^), (4.28a) 

ai2(k,u;) = 2^ 21.2 (4-28b) 

where we have restored the canonical dimensions and 
reintroduced the mass m. The leading order of the ex- 
pressions (4.28) (i.e., the first term) coincides with that 
found by Gavoret and Nozieres.^ The additional terms 
proportional to Qn give the next-to-leading contribution 
to the Green functions coming from the divergent dia- 
grams. These diverging sub-leading contributions were 
identified for the first time in Ref. 8 without calculating 
the coefficients. Later, by improving the Bogoliubov ap- 
proximation and exploiting the analogy with spin-wave 
theory, Weichman found an explicit expression for Qu 
valid for the weak-coupling Bose gas.^^ Assuming a Lan- 
dau quantum- hydrodynamic Hamiltonian, the authors of 
Ref. 24 have derived the same expression given by (4.27) 
for the longitudinal correlation function. 

D. Extension to all orders of PT 

Although the IR behavior of Qjj has been found only 
at the one-loop order, we will show now that it is actually 
correct to all orders. We begin by studying the nature 
of the fixed point. The scaling of the running couplings 
near the fixed point is given by Eqs. (4.23) and (4.22) for 
vii and by the WI for vm and vutt'- 

vittip) = vu{p)/^io{p) ~ vii{p) (4.29) 
vtttt{p) =3vitt{p)/^IJio{p) ~ vii{p)p-^ . (4.30) 

One thus finds that vut always vanishes near the fixed 
point. For d = 3 Vutt flows also to zero, while for d < 3 it 
flows to a finite fixed point value v^^^^ . We can thus draw 
one important conclusion, namely that for e = the the- 
ory is asymptotically free. Marginal interactions become 
marginally irrelevant (they vanish like 1/lnp) and the 
results obtained at one loop are exact, in agreement with 
Benfatto calculations.^^ 

Unfortunately, for e > the theory is not free since 
Vtttt fiows towards a finite value. This fact may, in princi- 
ple, invalidate the results found at one-loop level, because 
higher-order terms are not negligible near the fixed point. 
Note, however, that contrary to what happens in critical 
phenomena, the asymptotic behavior of vu does not de- 
pend on the position of the fixed point (w^t^t) but only on 
its very existence. This can be verified from Eq. (4.30). 
When Vttttip 0) is finite, vii{p) ~ and the exponent 
does not depend on the values of interactions. For this 
reason it is possible to establish the scaling of vu without 
a precise knowledge of the corresponding KG equation. 



Actually its "form" is enough, as we shall argue and ex- 
plain in the following by analyzing the PT. One has to 
prove either that vii{p) ^ p^,or the existence of the fixed 
point, since the two statements are equivalent. We will 
prove the former statement in the following. 

Let us consider the following two-loops diagrams con- 
tributing to F;;: 

-<i>- -0--00- 

These diagrams, once expanded in e, develop 1/e poles, 

and thus they contribute to the equation for vu [p) . With- 
out explicit calculations, we can infer the form of this 
contributions. We have seen that for p ^ Q, Gtt(k) and 
Git{k) do not depend on the explicit scaling of vu, thus 
they will not contribute to (4.31) with factors vu (at least 
at leading order as p ^ 0). This means, in particular, 
that the first diagram contributes to the beta-function 
only with a term proportional to vf^^, coming from the 
vertices. Its contribution is thus negligible with respect 
to the one-loop contribution (4.22). These classes of di- 
agrams do not affect the RG equations. 

Let us consider now the second diagram in (4.31), for 
which the situation is different since Gu is proportional 
to l/vii{p), but each internal Gii must appear in combi- 
nation with two vitt- From the scaling form at one-loop 
level, that we assume valid for the moment, we find that 
as p — > 0: 

-f«(p)e«(fc)--f#--%- (4.32) 
vii[p) 3 

Taking into account Eq. (4.32) this implies that the form 
the contribution of the second diagram in (4.31) is of the 
same form of (4.22), i.e., it has exactly the same p depen- 
dence. It does thus not change the form of the solution 
for Vu, even though it can give a multiplicative factor in 
the final scaling. 

The last diagram in (4.31) has a behavior similar to 
the second one. It depends on vu only through the two 
vitt connected to the external Gii- The internal part gives 
a finite contribution in p and again, once added to (4.22), 
it does not change the scaling of vu. We thus conclude 
that the contributions of aZZ two-loops diagrams to (4.22) 
do not change the IR behavior of vu. 

The same arguments can be extended at higher order 
in loops. The structure of the PT remains the same. The 
presence of even only one Git is enough to suppress the di- 
agram due to the vm non compensated by any l/vu. The 
Gil appears always in combination with two vut, giving a 
finite non-dangerous term. All diagrams thus are either 
negligible, or give a finite contribution for p ^ multi- 
plied by two Vitt coming from the two vertices connected 
to the external lines. 

The RG equation for vu to all orders has thus the fol- 
lowing simple form: 
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Provided the sum of the series does not change sign with 
respect to the sign of the first term, for p — > we obtain: 

vuip) = ^f^^ep^ (4.34) 

where a < 1 is the value of p for which the leading be- 
havior established by (4.33) holds and A{vff^f) is the sum 
of the series of Eq. (4.33). This concludes our proof. 
Higher-order diagrams do not change the IR behavior, 
as they only can affect the finite coefficients. 

We have thus proved that the scaling (4.24) holds ap- 
parently for any value of e > 0. One has to recall that 
the procedure breaks down at d = 1, since the conden- 
sate vanishes in that case, corresponding to e < 2. In 
conclusion, we have proved that the IR behavior found 
at one loop is maintained for 1 < d < 3. 



E. A non-trivial fixed point without anomalous 
dimensions 

We have shown that the fixed point is characterized 
by a finite value of the interacting constant among trans- 
verse fluctuations, similarly to what happens for the (j)^ 
theory. We have also shown that no anomalous dimen- 
sion appears in the transverse correlation function. This 
is clearly due to the constraints originating from the WI. 
The question is now how these constraints can be im- 
plemented in the structure of the perturbation theory. 
Actually, the same diagrams responsible for the anoma- 
lous dimension in the (j)'^ theory are also present here in 
Ttt- For instance: 

. (4.35) 

These diagrams, however, do not contribute in practice 
because other diagrams contribute with opposite sign and 
remove the divergence. This can be realized to all orders 
by looking at the form of Eq. (4.32). The structure Qiivf^^ 
acts as an effective vutt interaction mediated by the longi- 
tudinal field. According to Eq. (4.32), it has the opposite 
sign with respect to the fixed point interaction v^^^^. All 
contributions to the effective interaction are given by the 
following graphical equation: 

vtU = X + + ^ + ^ • (4-36) 

Evaluating the diagrams entering Eq. (4.36), one finds 

that at the fixed point v^l^l vanishes due to the exact 
cancellation of the fixed point v^^^f with the interaction 
mediated by Qu. 



This explains why no (/)'*-like divergences appear in the 
contribution to Tn- This does not mean, however, that 
the theory is free at the fixed point. As a matter of fact, 
Vf/fl is not defined in terms of one-particle irreducible 
diagrams for Tutt- Thus the above cancellation is only 
partial in some diagrams for Tu. This can be seen by 
expanding in a typical diagram for Tu the definition of 
v^/J. (represented here by the dashed blob): 

-co-=-oo-^-0- 

+0-0 ("'') 

The above diagram vanishes, but the last diagram does 
not contribute to Tu since it is not one-particle irre- 
ducible. Thus, vtttt does affect the equation for vu even 
if it gives only finite renormalization, as shown before. 
On the other side, one can verify that all vtttt terms 
appearing at all orders in Ttt are exactly canceled by 
the mediated interaction. We have thus explicitly shown 
how the diagrams responsible for the appearance of the 
anomalous dimension in 0^ theory get eliminated in the 
present case. This is the counterpart of what happens in 
the phase- amplitude variables, where the integration of 
the amplitude fields reconstructs an effective full-gauge- 
invariant action. 



F. IR behavior of the the response functions 

In Section IV C we have found the IR behavior of the 
quadratic Green's functions. Only Qu is affected by the 
nontrivial scaling of the four quadratic running couplings. 
We found that the form of the spectrum is unaffected by 
the RG flow. In the same way, we find that the phys- 
ical response functions, speciflcally the density-density, 
density- current, and current-current response functions 
{Gfi.v{k)) depend only on RG-invariant combinations of 
the running couplings. This is not surprising, since the 
expressions found for these functions in the limit k ^ 
by Gavoret and Nozicres depend only on two flnite phys- 
ical quantities, n and c^. As an example, we consider 
explicitly the density-density correlation function Q-flo- 
We obtain its expression from the general relation: 

Q-Mp) = r;.^(p) + ^jA-p)Sij{-p)ri-Ap) (4.38) 

with n = v = 0. The scaling of the composite F func- 
tions has been linked to the quadratic running couplings 
through the WI (3.26) and (3.28): 

Fi;o(k, u>) ~ witnl^"^ , Ft;o(k, w) ~ Uttul^^uj . 

Substituting this expression into Eq. (4.38) and using 
Utt —^-,00 /no in the IR, we recover the Gavoret- 
Nozieres result (with the mass reinserted explicitly): 
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nkV(nom) 

QMP) = 2 I 21 2 



(4.39) 



No anomalous terms appear. The same is true for the 
other response functions. 



V. CONCLUSIONS 

Infrared divergences plague perturbation theory of the 
interacting Bose gas without the presence of a true phys- 
ical instability for the system. This means that IR diver- 
gences must disappear from the final physical quantities, 
like the thermodynamical derivatives. However, it is not 
trivial to control the divergent parts of the perturbation 
theory whithin approximated shemes. Divergences are 
due to the strong fluctuations of the phase of the broken- 
gauge-symmetry order parameter. In this paper, we have 
set up a RG approach to treat these singularities in a 
systematic way. Since they must disappear in the final 
results, most divergences are actually related. We imple- 
mented these connections through Ward identities. We 
found that WI are sufficient to guarantee the stability of 
the superfluid phase with respect to phase fluctuations. 
As a consequence, the sound velocity of the macroscopic 
and microscopic excitations is an invariant of the RG 
flow. 

Nevertheless the divergences are not totally ineffective. 
The WI allow to reduce the set of independent running 
couplings to a single one, for which we have derived a 
flow equation at the one-loop order. Within this approx- 
imation, we found that the theory becomes free at d = 3, 
i.e., all running couplings flow to zero in agreement with 
Benfatto result. "'^^ One-loop order proves thus enough to 
obtain the exact result at d = 3, that gives a logarithmi- 
cally divergent longitudinal correlation function Qu. For 
1 < d < 3 the fixed point is characterized by a finite 
value of the coupling Vutt among transverse fluctuations, 
being a truly interacting fixed point. Even in this case we 
found that the one-loop result gives the correct exponent 
for Qii. Higher-order diagrams are not negligible like as 
in the = 3 case, but they contribute like the one-loop 
ones, and thus they do not change the IR behavior. Con- 
cerning the two other correlation functions Qu and Gtt, 
wc found that their IR behavior is determined by WI and 
no anomalous dimension appears. This result can be sur- 
prising for the transverse one, since at the fixed point the 
interaction is finite and an anomalous dimension could 
appear. We showed, however, how the structure of the 
PT adapts in a nontrivial way near the fixed point, in or- 
der to keep Vutt finite and fulfill WI. The transverse field 
is effectively free once one adds to Vtttt the interaction 
mediated by the longitudinal fields. 

In conclusion, our approach has enabled us to prove, 
within a coherent frame, the disappearance of the diver- 
gences in the physical quantities while keeping a firm con- 
trol on singularities during the procedure. The nontrivial 
evolution of the structure of the propagator from large 



momenta (Bogoliubov) to small momenta (fixed point) 
can be readily followed in terms of the flow of the run- 
ning couplings. This approach can be useful, in general, 
to study stable systems where singularities appears due 
to a broken continuous symmetry without true instabili- 
ties. 
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APPENDIX A: GAUGE INVARIANCE 
CONSTRAINTS FOR THE RUNNING 
COUPLINGS 

The introduction in Eq. (2.34) of the couplings asso- 
ciated with the monomials of the xi/t fields may lead to 
some difficulties in defining gauge invariance itself. Since 
the theory we arc constructing must satisfy the original 
gauge invariance defined in Eq. (3.1) in terms of the tp 
fields, we need to specify the the equivalent transforma- 
tion for the action (2.34) in terms of the new fields. In 
this way, wc can establish in a different way a connection 
among the bare couplings of the theory. 

To discuss this point we proceed by steps. We begin 
with the global part of the action and introduce the fol- 
lowing simplified notation for the couplings: 

^ = E ^i^n^Wt'i^) + ^li^i + «o + "t(V't + «*)• 

^-^ mini 

n,m 

(Al) 

In this expression, ai/t are two arbitrary real numbers. 
The global gauge invariance in terms of these fields reads: 



R(x) 



ijji (x) + ai 
tjjt{x) +at 



Oil 

at 



(A2) 



where R is defined in Eq. (3.3) and the external fields 
change as follows: 



a;(-^) 



R 



(A3) 



Note that the transformation (A2) depends on ai/^. This 
is due to the fact that the new action is no longer in- 
variant under rotations around the origin in the space 
{tphtpt), but rather around the point (— a;,— Q/). It is 
also clear the the original action is manifestly gauge in- 
variant, since it is a sum of gauge invariant terms (essen- 
tially IV'P")- By contrast, the action (Al) is invariant un- 
der the transformations (A2)-(A3) only for specific values 
of the set of couplings {vnm}- It is easy to find these con- 
straints by implementing the transformation (A2)-(A3) 
infinitesimally: 
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- at)Sx 



(A4) 



Imposing the invariance of S given by (Al) under the 
transformation (A4) gives the following set of equations 
(with at =0): 

-m Vn+l,m-l + n Vn-l,m+l + Oil Vn,m+l = . (A5) 

One can readily verify that these conditions coincide with 
the WI for ki = 0, the two procedures being equivalent. 

At this point it is easy to establish a link between the 
old couplings {n,v} and the new ones {vn,m}- By us- 
ing Eq. (A5) repeatedly, we find a relation among all the 
running couplings and reduce the whole set to the inde- 
pendent ones. This can be shown by noticing that all the 
running couplings can be derived by knowing only v„fi 
for all n, since the coupling Vn,m is related to Vn-i,m and 
Vn+i,m-2 W (A5). Thus, for any positive integer k we 
can draw a line in the n — m plane that starts at m = 
and n = k and end at m — 2k and n — 0. If the run- 
ning couplings, identified by the pair (n,m), which lye 
between the axes and this line are specified, then given 
Wfe+1,0 one can obtain all the couplings along the line 
Vn,-2n+2k for < 71 < /c + 1, thus extending the por- 
tion of the plane where the Vn^m are known. Proceeding 
in this way, all the running couplings can be obtained if 
Vkfl is known. Thus Vkfl can be chosen to be the indepen- 
dent variable. The usual action is recovered by choosing 
Vifi = —2ai{^—vaf), V2.0 = —2{fi — 3vaf), and all others 
Vn,o = (for at = 0). All other couplings can be found 
by the use of (A5). 

We consider now the local gauge invariance. The in- 
variance (3.1) holds for the original action. We need to 
clarify the precise meaning of requiring that the action 
written in terms of the new couplings Vn,m is invariant 
under (3.1). First, we face the problem of defining a /y, 
coupling, since there is no more trace of it in the new 
action (Al). The implementation of the invariance (A5) 
will thus involve a more complicated transformation of 
the fields, since no simple coupling couples to ji/jp. To 
find the proper transformation of the couplings, we need 
to consider the transformation under an infinitesimal lo- 
cal gauge transformation of the derivative term: 



S iwit Jdxijjt{x)dTi>i{x) 

-i^ jdx [i,^ + ^2 + 2aii^i] dr5x{T). (A6) 
This implies the following transformation (with at = 0): 



^2,0 - iwit drSxir) 

vifl vifi - iwitai dr6x{T) . 



(A7) 



One can readily recognize in (A6) the field ji/'P writ- 
ten in terms of the new variables. An equivalent ap- 
proach to define the local r-gauge invariant is then to 



introduce a new external source fj,extWit/2 that couples 
to IV'P (which will be set to zero at the end of the cal- 
culation). In term of fiexL, gauge invariance is defined as 
usual by Hext ^ Mecct + dr x(r). 

We consider now the gauge transformations involving 
spatial variations. In that case we have the following 
gradient terms in the action 



5 j dx\V^t\'^ + zii j dx\V^i 
where infinitesimal variation gives: 



(A8) 



2 Jdx{ztt M^i'tjySxix) - zu i>tiyi^i)V5x{x)} 

+2{ztt-zu) Jdx{\/iJtWi)Sxix). (A9) 

The global gauge invariance of the action [with x{x) = x] 
gives the condition zu = ztt- In this case, the resulting 
variation of the action reads: 

2 Ztt j dx {(V'; VV't) - (V't VV'i)} Vx(^) (AlO) 

that corresponds to the field coupled linearly to A(x) in 
Eq. (2.1). 

In this way, we have clarified the precise meaning of 

local gauge invariance in terms of the new couplings. 
These results can be used, in particular, to determine 
the composite-field insertions of Appendix B. 



APPENDIX B: WARD IDENTITIES AT 
ONE-LOOP LEVEL 

In this Appendix, we compute explicitly the one-loop 
contribution of relevant vertex functions. We can thus 
analyze the validity of the WI through an explicit ex- 
ample. In particular, the cancellations of the divergent 
contribution in the WI will be demonstrated. 

We begin with the analysis of the WI (3.12) when the 
coupling vi is nonvanishing. We will start by considering 
all running couplings present in the original action (2.1). 
The one-loop diagrams contributing to F; and F^t are: 



....0 



4 



^vittJ2^tt{q) 



-^viii^gil{q) , 

1 

-^Vtttt'^Gttiq) : 



: Vim 




-vfttJ2^ll{k-q)Qtt{q) 
-vftt^Gtiik - q)Qit{q) ■ 



(Bl) 
(B2) 
(B3) 
(B4) 
(B5) 
(B6) 
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For non- vanishing vi the conditions on the bare running 
coupUngs can be obtained from the WI (3.12), (3.13), 
and (3.9) (once Fourier transformed). They give, respec- 
tively: 



■tploVtt -vi=0, 

vitt-ipio + vtt - vii = , 
vim^Ho + '2vitt - Via = , 



(B7) 
(B8) 
(B9) 



while Eq. (3.19) is not modified (cf. Eq. (3.15)). The 
validity of Eq. (3.12) at the one-loop order can be now 
verified. Since Fj is the sum of u;, (Bl), and (B2), 
while Ftt is the sum of vtt and (B3)-(B6), we can readily 
evaluate Eq. (3.12) for vanishing external momentum. 
The lowest order is clearly satisfied by the bare values 
vi = — 2a(/z — va-^) and Vtt = — 2(/U — va^). To verify 
that also the next order (one-loop) is correct, we can use 
the following identity 



Qli{q)Qtt{q) - Gu{q)Gti(,q) 



Gu{q)-Gtt{q) 
Vtt - vii 



(BIO) 



holding for utt = and ztt = zn. Adding all the contri- 
butions and using (BIO), one can verify that the WI hold 
exactly. 

However, in obtaining a renormalizable theory, all the 
irrelevant running couplings have been discarded (for in- 
stance, viii or viitt)- As anticipated in Section III A, the 
WI are invalidated by this fact but they continue to hold 
at the leading order for A; ^ 0. This can be verified 
by considering the contributions to F; and Ttt when the 
marginal running couplings are set to zero. In this case, 
the contributions are (B1),(B3), (B5), and (B6), and the 
WI is not satisfied. We cannot test directly on this WI 
the cancellation of the leading divergent terms, since F; 
and Ftt(O) are finite. One could, however, verify that the 
ultraviolet cutoff divergence is still perfectly balanced in 
the WI; more instructive examples can be obtained by 
looking at higher-order WI. 

Let us consider accordingly the one-loop contributions 
to Til, ^itt, and Ftttt due to the marginal couplings (from 
now on we set to zero all the irrelevant running cou- 
plings). On the same footing, the Green's functions are 
given by the expression (2.36) with zu = 0, and Utt in gen- 
eral nonvanishing. The mean-field condition is applied to 
the lowest-order values for vu and Vtt, which thus vanish. 
Notice that, using the above conditions, (BIO) becomes: 



Gtt{q)Gii{q) - Git{q)Gti{q) = - 



Gtt{q) 



vii 



(Bll) 



The one-loop diagram contributing to F;; (fc) gives: 

--Cy- =-'^T.^tt{k-q)Gtt{q). (B12) 

We evaluate also the one-loop contributions to 
ritt(fci) ^2, ^s)- For simplicity, we write only the value 



of the diagrams for vanishing external momenta, since 
the value for finite k can be reconstructed from the fol- 
lowing expressions by inserting the correct momenta in 
each internal Green's functions and adding the contribu- 
tion of the permutations of the external arguments. We 
report below all the one- loop diagrams for Tut'. 



(B13) 



2 - 



vltY.^tl{qfGtt{q) , 
^4tY.^tt{q)Gtl{q), (B14) 

Q 

-vfttJ2^tt{qfGu{q) (B15) 



Vitt^tttt 



Y^Gttiqf. (B16) 



Only in diagram (B14) there are two different permu- 
tations. With the same notation, we consider now the 
one-loop contributions to Ftt^j. We write the following 
equation: 



Ttttt = 6 



+3 



+ 12 



+ 6 



S. 



(B17) 



The numbers that multiply the diagrams indicate the 
non-equivalent permutations of the external ki . The con- 
tribution of these diagrams is identical to the one given 
in (B13)-(B16), apart from the number of permutations 

and an overall factor vuti/vm to be multiplied. Instead 
S represents the sum of the following set of diagrams: 



24 ! ' 

/ — S 
/ \ 



= -^vt,,Y.Gu{q)' 



(B18) 
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= -l2vtttY,Git{qfGn{q)Gtt{q), (B19) 



= ^vi,Y.Qit{qfGu{q)Gtt{q), (B20) 



-^vtuY.Gtt{qfGu{qf. (B21) 
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Adding all these contributions together and using (Bll), 
we finally obtain for the F functions at vanishing external 
k: 



^tttt — ~ 



(B22) 
(B23) 



The above expressions are divergent when d < 3. An IR 
cutoff should thus be introduced, with the vertex func- 
tions diverging when this cutoff goes to zero. It is read- 
ily verified that the leading diverging contributions of 
each vertex function simplify exactly when (B22), (B23), 
(B12), and (B3) are inserted into (3.13) and (3.15). Note, 
however, that there exist finite terms (specifically, the 
contribution from Ttt) that do not vanish, so that the 
WI is not exactly satisfied. This fact had to be expected 
and docs not spoil the RG identity among the running 
couplings, since within the scheme of the dimensional 
renormalization they are valid for the diverging contri- 
butions. The PT is thus preserving the gauge invariance 
at the leading order in the IR divergences. 

We consider now Eq. (3.11) as the simplest example of 
cancellation of the diverging coefficients of the external 
frequency kg- To verify this, we need to evaluate Tu and 
r;;o for finite kg- The one- loop contribution is: 



Titik) 



-vitt 



^Gtt{k- q)git{q) , 



(B24) 

ri;o(fc) = - - -^vittwit'^gttik - q)gtt{q) ■ 

(B25) 

When we evaluate Eq. (3.11) for ki = and finite ko, 
wc obtain (we have used Qit{k) = Qtt{k)kQWit / vu from 
(2.36)): 

-\vittWuY,Qtt{k - q)Qtt{q){qo - ko/2) . (B26) 



In order to extract the /cg coefficient, let us perform the 
change of variables q = q' -\- k/2. Since an ultraviolet 
cutoff A is assumed, the above change of variable should 
change the limits of integrations. We neglect this fact, 
since the difference is of higher order in go/A, and we arc 
interested in the go ^ limit. With the above change of 
variable, we readibly obtain that the leading contribution 
of (B26) vanishes. 

It is not difficult to evaluate also Ti.^i{ko,yi) and insert 
the result into Eq. (3.11). We can evaluate then all T 
functions at finite values of ko and ki. The above results 
remain valid for the first part of the WI, since no use was 
made of the fact that kj was vanishing. However, the 
contribution of the vector part of the WI is sub-leading 



with respect the one calculated in (B26). To check the 
vanishing of the vector part, one should thus perform 
the above calculation by keeping also the next-to-leading 
order. 



APPENDIX C: ONE-LOOP RG EQUATIONS 

In the text we have shown that a single coupling is 
actually scaling independently. In this Appendix, we ob- 
tain the one-loop RG equations for the four marginal 
couplings entering the propagators {vu, wn, uu, and zu) 
and verify that the exact constraints found in Section 
IV B are fulfilled. We need to calculate explicitly the one 
loop contributions to Tu [depicted in (B12)], Tu [depicted 
in (B5)], and [depicted in (B5), and (B6)]. Accord- 
ing to the general procedure established in Section II B 
no correction to a need to be considered. The integrals 
entering these expressions can be written in terms of the 
following integrals (when the expression for the propaga- 
tors in the small- momentum limit is considered): 



^n,m (^) 



d'^+^q g;^(q2)'"/2 



{2TrY+^ {k + qfq^ ' 
We thus have at the one loop order: 



Vu{k)=vn- 



/o,o(fc) 



tt 



(CI) 

(C2) 
(G3) 



r;t(A;) = witko + ^"^2 -^i,o(fe) 
Vtt{k) = uttkl + ztt^i^ - 2 [i'^tmi + wl)l2,o{k) 



-ztmihA^) + wf^knIlfl{k)\ 



(G4) 



The In,m integrals are performed by standard 
procedures: 



dx 



r(l)2 J (27r)'^+i Jo [x{k + qy + [l - x)q^Y 



- I dx 



d'^+'^q f{q- kx) + f{-q- kx) 
2io"",/ (2^0^ [g2 + k^x{l - x)]^ 

where f{q) = g"q"*. One can now use the identity 

T{{d+l)/T))T(a-id+l)/2) 



+°° x'^dx 

(1+X2)« 



2r(a) 



(G5) 



to obtain the; analytic expression to be continued to com- 
plex values of the variable d. As a matter of fact, the 
integral (C5) converges only for d — 2a > 0, but its ana- 
lytic continuation given by the Eulcr F-functions on the 
right-hand side of (C5) is defined for any value of the 
complex variable d, except for a discrete set of points 
where simple poles arc present. Dimensional rcgulariza- 
tion proceeds by subtracting these poles that in our case 
are in 1/(3 — d) = 1/e. We thus expand the integrals in e 
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keeping only the leading order. This gives the following 
expressions for Fij (we have absorbed the phase-space 
factor l/(87r^) in the definition of the running couplings 
m vii/Sn'^ and wu wu/VStt^)-- 



^u{k) = vii 



4 ik'f' 



2 ^2 
tt 



witvu , (k ) / c 
ko 



<ki {ey" 



where c? indicates the combination of running couplings 
forming the square of the sound velocity [cf. Eq. (3.38)]. 
As a first result, we note that no divergence appears in 
the dependence. This implies that zu does not scale, 
as demonstrated using the WI in Section IIIB. The RG 
equations obtained by minimal subtraction for the re- 
maining three couplings read: 



^dm ^ vfiip) cjp) 



dp 

dutt 
dp 



2V'fo(p)4 Co 




vii{p)wu{p) f 

2V'fo(p)4 V 


c_{py 

Co . 


Kip) ( 


'<p)_ 




. Co 



(C6) 
(C7) 
(C8) 



where ipio{p) = ipio{^)p'^^'^ due to its bare dimension. It 
is now possible to veriiy that the RG equations admit the 
invariants witip)/vii{p) and c?{p), by simply performing 
the derivative with respect to p and substituting the RG 
equations. This shows at the one-loop level the validity of 
the identities proved at all orders in Section IV B. These 
identities can be used to eliminate two out of the three 
running couplings. For the last one (w;;, for instance), we 
need to solve the RG equation explicitly. The solution is 
given and discussed in Section IV C. 
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